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In a previous work [Q] , we have shown that a spinor condensate confined in a periodic or double- 
well potential exhibits ferromagnetic behavior due to the magnetic dipole-dipole interactions be- 
tween different wells, and in the absence of external magnetic field, the ground state has a two-fold 
degeneracy. In this work, we demonstrate the possibility of observing macroscopic quantum spin 
tunneling between these two degenerate states and show how the tunneling rate critically depends 
on the strength of the transverse field. 
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Tunneling, a process in which a system penetrates into 
a classically forbidden region (e.g., a potential barrier), 
is an intrinsically quantum effect with no classical coun- 
terpart. Of all tunneling effects, macroscopic quantum 
tunneling (MQT), the tunneling of a macroscopic vari- 
able of a macroscopic system represents a particu- 
larly intriguing and interesting scenario as it touches the 
boundary between the classical and the quantum world 
and may help shed light on the quantum-classical inter- 
face. As such, it is one of the most striking manifestations 
of quantum mechanics. The tunneling of large magnetic 
spins has recently received much attention, both theo- 
retically and experimentally Q in view of its promise 
as one of the few realistic candidates for an experimen- 
tal demonstration of MQT, and also because of its con- 
nection to quantum computing. Despite considerable ef- 
forts, though, there are not many clear and definitive 
demonstrations of MQT in spin systems so far — al- 
though there have been some indications that tunneling 
might be the underlying reason for some observed re- 
sults 1^, 1^. The reasons for this difficulty are as follows: 
First, most theories apply to single system, while the 
conventional magnetic materials used in the experiments 
contain many domains, each possessing its own set of pa- 
rameters such as magnetic anisotropy and barrier energy; 
second, due to the difSculty of cooling the samples down 
to ultracold temperatures, thermal processes cannot be 
completely excluded; and third, the spins in solid ma- 
terials are inevitably imbedded in a crystal matrix and 
the spin-matrix interaction |^ complicates the physical 
picture. For these reasons, macroscopic spin tunneling 
remains one of the most anticipated, yet elusive, quan- 
tum phenomena. 

In this paper, we demonstrate the possibility of observ- 
ing spin tunneling in a spinor atomic condensate trapped 
in a double-well potential, thereby eliminating most of 
these difficulties. We remark at the outset that while the 
inter-well tunneling of condensates has been previously 
considered, all previous studies focused instead on the 
tunneling of an external degree of freedom, the conden- 
sate center-of-mass motion [[zf p|, p|, p^. 



In an earlier paper , we showed that because of the 
long range magnetic dipole-dipole interaction, a spinor 
atomic condensate trapped in a one-dimensional periodic 
lattice potential or a double-well potential behaves as a 
ferromagnet. In the absence of external magnetic fields, 
the ground state spins of the "mini-condensates" con- 
fined in individual wells all align parallel to each other 
and along the lattice direction (the z-axis), giving rise 
to a spontaneous magnetization along z. For the sake of 
simplicity, we restrict the present discussion to a double- 
well potential, each well containing N spin-1 condensate 
atoms. In the tight-binding approximation [Q, they are 
described by the zero-temperature spin Haniiltonian ||l[ 

H = \'^{?,l + Sl)- 2,\SIS'^ + ASi • S2 - h{Sl + 5^), (1) 

where is the total spin of the condensate in the i^^ well, 
and S^'^'^ its cartesian components. The first term in 
(|l|) represents the on-site Hamiltonian of the spinor con- 
densate. It includes short-range nonlinear spin-exchange 
interactions , where the parameter X'^ is related to the 
s-wave scattering lengths |l2l and needs to be negative. 
The second and third terms in (|l|) arise from the site-to- 
site dipole-dipole interaction where A = ^g^Q/{Anr^) 
for pure magnetic dipolar interaction, with jb being the 
gyromagnetic ratio, /iq the vacuum permeability and r 
the distance between the two wells. The value of A can 
be greatly enhanced by the light-induced optical dipolar 
interaction if one chooses appropriate laser fields to form 
the potential well The last term describes the 

effect of an external transverse magnetic field, taken to 
be along the x-axis without the loss of generality. Here 
h ~ IbB, with B being the strength of the applied field. 
In the following, we assume that the nonlinear short- 
range atom-atom interaction is strong enough that the 
first term in Hamiltonia n (p] ) dominates over the mag- 
netic dipolar interaction | [15[ . As a result, the total spin 
quantum number for each mini-condensate is fixed to its 
maximum value N — the number of particles in each 
well. We can therefore neglect the first term in the Hamil- 
tonian (|^), since it is a constant of motion and commutes 
with the remaining terms. 
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FIG. 1: Classical ground state spin orientation, (a) For < 
h < 3NX, the ground state has a two-fold degeneracy, (b) 
For h > 3A''A, the degeneracy is removed and the spins are 
polarized along the transverse field. 



Before discussing quantum mechanical spin tunneling, 
let us first investigate the classical situation. The Hamil- 
tonian in that limit is still given by (|l]), except that the 
spins are now c-numbers that can be represented by vec- 
tors of fixed length N in spin space. For zero field, h = 0, 
it is easy to see that the classical ground state is two- 
fold degenerate with S = Si = S2 pointing along ei- 
ther the z- or (— z)-axis. Under the influence of a weak 
transverse field along the a;-direction, the two ground 
states move away from the z-axis and towards the x-axis 
while remaining in the a;2;-plane, as shown in Fig. For 



0<h<h. 



3NX, the two minima are located at 



= 7r/2±cos"i (h/hi 



where 9 is the angle between S and the z-axis. For h > 
hi''\ the two minima merge along the x-axis, and the 

(c) 

degeneracy is removed. Hence he can be regarded as the 
classical critical field strength, beyond which the system 
is completely polarized by the external field. 

Let us now turn to a quantum mechanical description 
of the system. Our goal is to investigate whether or not 
tunneling is present in the classically degenerate regime 
(i.e., when < h < hi'^^). We will present both a full 
numerical calculation and analytical results using the in- 
stanton technique. 

A well-known consequence of the tunneling between 
two degenerate states is the lifting of their degeneracy: 
The two new eigenstates are a symmetric and an anti- 
symmetric superposition of the original states character- 
ized by an energy difference (or tunneling splitting) Ae 
inversely proportional to the tunneling rate. The quan- 
tity of interest to determine the occurrence of tunneling 
is therefore this energy difference between the two lowest 
eigenstates of the Hamiltonian. We determine it by ex- 
panding the Hamiltonian (^ onto the basis spanned by 
S'f (8) S'l, and evaluate numerically the eigenvalues of the 
resultant {2N +1)^ x {2N + 1)"^ matrix. Fig.|(a) summa- 
rizes the result of this analysis. It shows that Ae is essen- 
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FIG. 2: (a) Energy splitting Ae as a function of the transverse 
field strength. The symbols represent numerical results and 
the solid curves are analytical results obtained from Eq. 
(b) The quantum critical field strength hi''^ as a function of 
A^. hi''^ is defined as the value of h at which Ae — 0.1 A. 



tially zero for small values of h, but becomes finite when 
h exceeds a threshold value hi"^^ — the quantum critical 
field strength. This means that for < h < h^c \ quan- 
tum mechanics agrees with classical mechanics in that 
the system is degenerate. However, for h'^c'^ < h < hi^\ 
even though the system is still degenerate in the classical 
picture, the presence of tunneling removes the degener- 
acy in the quantum treatment. Fig. ^(b) displays h'f ^ as 
a function of N, from which one sees that h'f ^ increases 

(c) 

with N and approaches he as N tends to infinity. In 
other words, as TV increases, the system, as expected, 
behaves more and more classically. 

To gain some analytical insight, we first notice that 
Hamiltonian dll) can be rewritten as 



-S2 - hS^ 
2 



where S = Si S2 and S' = Si - S2. For h < h^c \ Si 
and S2 are tightly bound together, such that S'^ ~ and 
S is approximately a constant of motion with quantum 
number S = 2N. Hence, neglecting the constant terms, 
the effective Hamiltonian of the system reads 



^ -^{S'f ~ hS^ 



(2) 



Hamiltonian (^ describes a quantum spin with the easy- 
axis anisotropy in a transverse field and it has been ex- 
tensively studied in the context of spin tunneling . Us- 
ing the instanton technique, which has been proved to 
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be quite accurate for S* > 5, the tunneling splitting be- 
tween the two classically degenerated ground states can 
be expressed as 

Ae=pc^V5,/(27r)e-^% (3) 

where p is a prefactor on the order of unity [by fitting 
the numerical results with Eq. (ph, we find p « 2.75], 

(c) . . . 

uj = he X is the typical instanton frequency with x ~ 

\Jl- {h/h\f^Y, andS'c = 27Vln[(l + x)(l-a;)]-4iVa; is 
the classical action. The solid curves in Fig. |^(a) repre- 
sent the values calculated using Eq. (|), which agree well 
with the numerical results. 



A more intuitive understanding of this threshold be- 
havior in quantum tunneling can be obtained from a 
quantum perturbation viewpoint. Tunneling results from 
the transverse magnetic field. If we regard the last term 
in Hamiltonian (|l|) as a perturbation to the rest of the 
Hamiltonian, which possesses two degenerated ground 
states \N^N) and | — N,—N), where we have labelled 
the states with the eigenvalues of Sf and S"!, the tun- 
neling level splitting is then, in the high order perturba- 
tion theory given by the shortest chain of matrix 
elements and energy denominators connecting these two 
states For the situation at hand, this chain can be 
represented as 



\N,N) ^ -^{\N-l,N) + \N,N^l))^\N-l,N-l)^^{\N-2,N-l) + \N-l,N-2)) 
\/2 V 2 

...^ -^{\~N,^N+l) + \-N + l,-N))^ \-N,-N). 
v2 



r 



Since the spin has to travel a distance of iN to reach 
from one ground state to the other, the tunneling appears 
minimally in the (4A^)*^ order of perturbation theory. 
The energy splitting associated with this chain can be 
calculated as 



Ae = 16A 



8AiV' 



2V2A 



4JV 



r (3/2 + + 1/4 - 
r (iv + ^iV2 + 1/4 - 1/2) 



\aV2nxJ 



where the last line is obtained under the limit N 1. 
The power law dependence on the field strength indicates 
that Ae is vanishingly small for large N unless the term 
inside the bracket exceeds one, resulting in a critical field 
strength at (4\/2/e)A^A « 0.7hi''\ Ahhough this treat- 
ment explains the threshold behavior, it fails to predict 
accurately the critical field strength. In particular, it 
fails to predict the A^-dependence of the ratio h^c^ /h\^'' 
see Fig. 1(b)] . This failure can be easily understood 
since for h close to the critical point, it is no longer valid 
to regard the transverse field as a perturbation. 

We now turn to the detection of tunneling. It is not 
practical to directly measure the energy difference be- 
tween the two lowest energy states, since Ae is relatively 
small compared to the total ground state energy. Instead, 
we use a method which is perfectly suited for the situ- 
ation at hand, where thermal effects are negligible, but 
may not be appropriate for more traditional solid-state 
systems. Our proposed detection scheme starts with the 
system prepared in one of the degenerated ground state. 



say, |A^, A^). The external field h is then slowly ramped 
up from zero to a final value, hj < he ■ Classically, the 
macroscopic spin would simply adiabatically follow the 
instantaneous ground state closest to the initial state. 
Quantum mechanically, though, this is only true as long 



as hf < hi'^K But for hf > he'", the system is in the 
quantum tunneling regime and does not reach an equi- 
librium. Rather, it oscillates back and forth between the 
two classically degenerate states, a phenomenon termed 
as macroscopic quantum coherence. The observation of 
such oscillation will be a direct signature of quantum tun- 
neling, but it requires the effect of dissipation to be small. 
Otherwise, spin relaxation will result, which is normally 
the case for the experiment on solid magnetic materials. 
The long coherence time associated with the atomic con- 
densate, however, makes it ideal for this purpose. In fact, 
macroscopic center-of-mass oscillation has already been 
observed in ultracold atoms trapped in optical lattices [|| . 
Fig. ^ depicts the numerical simulation of the measure- 
ment scheme in both the classical and quantum treat- 
ments. The quantum mechanical results are obtained 
by time evolving the Schrodinger equation in spin space 
with Hamiltonian (|l|) , using the Crank-Nicolson method, 
while the classical results are obtained by solving the dy- 
namical equation 



(9) 



dS, ^ SH 



(4) 



where H is the energy functional and has the same form 
as Hamiltonian but with the spin operators treated 
as c- numbers. Fig. ||(a) is for /i/ < h'fK There is no 
qualitative difference between the quantum and classical 
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FIG. 3: Time evolution of the expectation value of . Here 
N = 10. The solid lines represent the classical results and the 
dashed lines the quantum mechanical results. The strength 
of the transverse field, h, is ramped linearly from to /i/ 
from t = to t = 0.5, then stays ai hf afterwards. For (a), 
hf = 0.9hc, no tunneling is present; for (b), hf — 1.15/ic, 
macroscopic quantum coherence due to spin tunneling can be 
observed in the quantum mechanical results. The horizontal 
fines represent = ±2N^1 - [/i//(3iVA)]2, values of in 
the classical ground state with h = hf. The units for time is 
h/X, which is on the order of a few seconds for typical atomic 
parameters for pure magnetic dipolar interaction. 

results, a steady state is reached in the end (small os- 
cillations around this steady state persists due to nona- 
diabiticity). Fig. ^(b) is for hf > hi'^K The classical 
result is quite similar to that of Fig. ||(a). However, the 
quantum calculation clearly shows the oscillations of the 
system between the two macroscopically distinct states. 
These two states differ by a minus sign in the expectation 
value of , which can be easily measured experimentally 
using, for example, Stern-Gerlach technique 

In conclusion, we have shown that the quantum macro- 
scopic tunneling of spin is possible in a spinor condensate 
trapped in double-well potential. Compared to the more 
conventional solid state magnetic materials, our system 
possesses several decisive advantages. First and perhaps 
foremost, this is an exceedingly clean system character- 
ized by a few simple parameters, without the compli- 
cation of domain separations, and with well-understood 
microscopic physics; It is amenable to exquisite experi- 
mental control; A single parameter — the transverse field 
strength — is capable of switching on and off the tun- 
neling. Typical temperatures of the atomic condensate 
is in the nanoKelvin regime, hence the thermal activa- 
tion normally presented in the solid materials can be 
safely neglected Due to these reasons, we believe 

that this system is indeed ideal for studying quantum 



magnetism in general |T^, and macrosopic quantum 
tunneling in particular. Finally we remark that although 
we have considered a double-well potential in this work, 
we expect similar behavior for a condensate confined in 
a one-dimensional periodic lattice potential. 
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